Di-electron production from vector mesons with medium modifications in heavy ion 

collisions 
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We reproduce the di-electron spectra in the low and intermediate mass regions in most central 
Au+Au collisions by the STAR Collaboration incorporation of the STAR detector acceptance. We 
also compare our results with PHENIX data constrained by the PHENIX acceptance. We include 
the medium modifications of vector mesons from scatterings of vector mesons by mesons and baryons 
in the thermal medium. The freezeout contributions from vector mesons are also taken into account. 
The space-time evolution is described by a 2+1 dimensional ideal hydrodynamic model. The back- 
grounds from semi-leptonic decays of charm hadrons are simulated by the PYTHIA event generator 
and corrected by the nuclear modification factor of electrons from charm decays. It is difficult 
to extract the thermal contributions from those from charm decays in the invariant mass spectra 
alone and in the current detector acceptances. Other observables such as transverse momenta and 
collective flows may provide additional tools to tag these sources. 



I. INTRODUCTION 

The electromagnetic probes such as photons and dilep- 
tons are expected to provide clean signatures for the 
quark gluon plasma (QGP) in heavy ion collisions due to 
their instant emissions once produced [l|-0j- These ther- 
mal photons and dileptons contain undistorted informa- 
tion about the space-time trace of the new state of matter 
formed in such collisions. The invariant mass spectrum 
is usually divided into the low, intermediate and high 
mass regions (LMR, IMR and HMR), based on the no- 
tion that each region is dominated by different sources of 
dileptons. In the LMR, M <1 GeV, dileptons are mainly 
from vector meson decays and may be related to chiral 
symmetry restoration [8Hl3j. In the HMR, M > 3 GeV, 
dileptons are dominated by the Drell-Yan process and 
quarkonium decays. In the IMR, 1 < M < 3 GeV, it 
was argued that dileptons from semi-leptonic decays of 
correlated open charm hadrons are dominant [l4j . 

The medium modifications of the p meson spectral 
functions are successful in describing the di-muon en- 
hancement in the LMR of the NA60 experiment at the 
SPS energy [H El El- The PHENIX and STAR col- 
laborations also observed such an enhancement in the di- 
electron spectra at the RHIC energy 

E1E3- T^ ther- 
mal quark-antiquark annihilation in the QGP phase is ex- 
pected to give a measurable signal in the IMR for the de- 
confinement phase transition at RHIC energy (l^. How- 
ever, in this mass region, the di-lepton yield from semi- 
leptonic decays of open charm mesons increases rapidly 
with the collisional energy. The single leptons from open 
charm mesons and their dynamic correlations are ex- 
pected to undergo medium modifications. The question 
is: to what extent the di-leptons from charm hadrons 
with medium modifications mix up with the thermal con- 
tributions from the QGP in the IMR. Another issue is 
that the dilepton spectra measured by the STAR and 
PHENIX collaborations are very different in the LMR. 
It is worthwhile to to look at this disagreement closely 



by using the Monte Carlo simulation incorporating the 
different acceptances of STAR and PHENIX detectors. 

In this paper, we try to reproduce the data of di- 
electron invariant mass spectra in the LMR and IMR in 
central Au+Au collisions at 200 GeV. We will include the 
medium modifications of the vector mesons and charm 
hadrons. The acceptances of STAR and PHENIX de- 
tectors are incorporated in our calculation. We will use 
a 2+1 dimension ideal hydrodynamic model to give the 
space-time evolution of the fireball, where the parameters 
are determined by fitting the data of transverse momenta 
of long life hadrons (pions, Kaons and protons). The 
spectra of charm hadrons (D°,D ± , D s and A c ) are given 
by a simulation of the PYTHIA event generator. We ne- 
glect the Dalitz decay channel for 7r°: tt° —> e + e~7 but 
include those for r\ and u>: rj 



D+ e 7 and u) — > 



The contribution from pion's Dalitz decay is mainly be- 
low and irrelevant to our current range of the invari- 
ant mass. 



II. PARAMETERS IN HYDRODYNAMIC 
MODEL 



We use a 2+1 dimensional ideal hydrodynamic model 
[T§ | to give the space-time evolution of the medium cre- 
ated in heavy ion collisions. We choose two types of the 
Equation of state (EOS) (H^H, S95P-CE (CE) with 
complete chemical equilibrium to very low temperatures 
and a wide range of phase transition temperatures from 
184MeV to 220MeV, and S95P-PCE (PCE) with partial 
chemical equilibrium below chemical freezeout temper- 
ature Tchem = 165 MeV (H. After kinetic freeze-out, 
we use the Cooper-Frye formula [H [24| to obtain the 
momentum spectra for each hadron species 
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set 


EoS 


T chem [MeV] 


Tf [MeV] 


e/ [GeV/fm 3 ] 


SI 


S95P-CE 




136 


0.12 


S2 


S95P-PCE 


165 


136 


0.275 


S3 


S95P-PCE 


165 


106 


0.12 



Table I: Parameter sets. We choose most central collisions 
with b = 2.4 fm. The initial conditions are chosen to be: 
r = 0.4 fm , e = 45 GeV/fm 3 (T =395 MeV). 



where E M denotes the normal vector of the freezeout hy- 
persurface, Tf is the kinetic freeze-out temperature on 
the freezeout hypersurface, rij is the phase space distri- 
bution function for the baryon/meson species i which can 
be Fermi-Dirac/Bose- Einstein distribution and gi is its 
degeneracy factor, and pt denotes the transverse mo- 
mentum. 

Beside the EOS, there are some free parameters which 
should be fixed in the hydrodynamical model, such as 
equilibration time To, the initial energy density eo (or 
initial temperature To) and the kinetic freeze-out energy 
density e/ (or Tf). At the RHIC energy ^s^n = 200 
GeV for Au+Au collisions [25|, (2(|, we constrain these 
parameters with STAR and PHENIX data for the rapid- 
ity densities of multiplicities, dNi/dy, and the px spec- 
tra for long-life hadrons (pions, koans and protons). The 
distribution of the initial energy density is determined by 
the Glauber model with 5% of the contribution from bi- 
nary collisions. We focus on most central collisions with 
impact parameter b = 2.4 fm throughout the paper and 
assume a system with vanishing net baryon number. 

To make a comparison, we choose the same initial con- 
ditions for two EOS (r = 0.4 fm , e = 45 GeV/fm 3 
(or Xo=395 MeV)), but freeze-out conditions are differ- 
ent due to different relations between the energy density 
and temperature. The same initial conditions imply that 
the entropy densities are chosen to be the same for these 
two EOS. We choose T f = 136 MeV (e } = 0.12 GeV/fm 3 ) 
for CE and roughly reproduce the pt spectra for the pi- 
ons and kaons, as shown in Fig. (pTa)). But the proton 
yield is under-estimated. For PCE, we choose the kinetic 
freeze-out occurs at either the same energy density or the 
same temperature as CE, i.e. (e/,T/) = (0.12 GeV/fm 3 , 
106 MeV) or (0.275 GeV/fm 3 , 136 MeV). The p T spectra 
with PCE are shown in Fig. ([Ub)). The parameter sets 
are listed in Tab. (JTJ) . 

The PCE scenario has a feature compared to the CE 
one pj|: dNi/dy and px spectra in the low pt range 
(pt < 1 GeV) are almost independent of Tf. However, 
the high Tf will reduce the dilepton production rate by 
shortening the evolution time. So we can regard Tf as 
a tuning parameter to the dilepton production rate from 
the in-medium vector meson decay. 

We give a few comments about the yields of vector 
mesons, whose multiplicity ratios, p/tt,uj/tt and </>/7r, are 
given in Tab. pi]). Most of the decays of the uj and <fi 
mesons take place after the kinetic freeze-out. The yields 
of uj and 4> are in good agreement with the data and will 
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Figure 1: (Color online) Transverse momentum spectra for n + 
(black-dashed), K + (red-dash-dot-dot-dotted) and p (green- 
dotted) for central collisions with centrality 0-5%. The data 
are from the PHENIX collaboration 0. (a) S95P-CE and 
(b) S95P-PCE EOS are used. In (b), the differences between 
T f = 106 MeV (thick lines) and T f =136 MeV (thin lines) in 
the low pt range are small. 



contribute to the dilepton rate. We see that the ratio p/n 
at the freeze-out is smaller than the data, but we note 
that most of the p mesons decay in the thermal medium 
earlier than the kinetic freeze-out whose contribution to 
the dilepton production dominates the dilepton spectra. 
Hereafter we will use the parameter set S3 of PCE as a 
default choice unless stated explictly. The comparison 
will be made with other parameter sets. 



III. DILEPTON EMISSIONS IN HEAVY ION 
COLLISIONS 



In the thermalized medium, hadron gas (HG) or quark 
2;luon plasma (QGP), the dilepton production rate per 
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sets 


7T + 


V 


p/TT 


oj/n 


<j>/lt 


SI 


285.7 


12.2 


7.91 x 1CT 2 


7.49 x 10~ 2 


1.86 x 10~ 2 


S2 


270.2 


23.6 


7.78 x 1CT 2 


9.85 x 10~ 2 


2.89 x 10~ 2 


S3 


261.7 


24.6 


5.50 x 10" 2 


9.82 x 10" 2 


3.13 x 10 -2 


PHENIX 


281.8 


18.4 


1.03 x 10 _1 


8.98 x 10~ 2 


2.14 x 10~ 2 


STAR 


327 


34.7 


1.69 x 10" 1 




2.65 x 10" 2 



Table II: Rapidity densities dNi/dy for meson/proton yields 
in most central collisions with b = 2.4 fm. The PHENIX data 
are taken from Ref. [l6|, |26| |. where they only have the 7r + 
data in most central collisions and their p data are from the 
fragmentation model. The STAR data are from Ref. [27l.[2§|. 
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Here mi is the lepton mass, a = e 2 /47r is the fine 
structure constant with the electric charge e for lep- 
tons, p = (poiP) = Pi + P2 is the dilepton 4-momentum 
and M = yjjfi, n B = 1/ (e p - u/T - l) (T and u are 
the local temperature and fluid velocity respectively) 
is the Bose distribution function, 11^, is the retarded 
photon polarization tensor from the quark or hadronic 
loop, and H R = ^iffi* . For the partonic phase, 
given by the Born term reflects the lowest order process 
qq — > 7* — > . For the hadronic phase, 11^ is further 



related to the retarded vector-meson propagator Dy with 



V = p, (J, (j> via Inin^ = — (e 2 m v /g v ) ImDy , where gv 
is the photon- vector-meson coupling constant in the vec- 
tor meson dominance model, and my is the vector-meson 
mass. The retarded vector meson propagator is 
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where 11^ is the contraction of the retarded vector meson 
polarization tensor. 

We use the hadronic many body effective theory (29| 
to calculate in-medium p spectral functions by scatter- 
ing with surrounding mesons. Though we assume a net 
baryon free system at RHIC energy, the effective chemi- 
cal potentials in the PCE scenario will give a considerable 
number of baryons (30j . and we assume that there are 
equal number of anti-baryons which give the same con- 
tribution as the baryons to the p spectral functions. To 
include the baryonic (including anti-baryonic) contribu- 
tions, we use the empirical scattering amplitude method 
(3l1 | which agrees with the hadronic many body effective 
theory (3^1 ■ Here we only consider the coupling of the p 
meson with baryonic resonances in the medium and set 
the momentum q = 300 MeV for the p meson in-medium 
propagator. The in-medium p meson spectral functions 
with and without baryonic contributions are shown in 
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Figure 2: The imaginary parts of the in- medium p meson 
propagators (or in-medium spectral functions) with (thick 
lines) and without (thin lines) baryonic contributions. The 
chemical potentials in PCE EOS are used. 



Fig. ([2]) at different temperatures but at a fixed mo- 
mentum q = 300 MeV. The differences between with and 
without baryonic contributions are larger at low temper- 
atures than at high temperatures. 

Since the collision rate in a meson gas around the tran- 
sition temperature indicates a large broadening of 4> me- 
son spectra due to binary collisions (33| . we include this 
effect via a schematic estimate as follows (with T = 150 
MeV), 
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In Fig. (H3a)) we show the invariant mass spectra of 
thermal di-electrons at the RHIC energy 200 GeV for 
most central Au+Au collisions. Beside the p compo- 
nent (red-dotted line) in the hadronic phase in an early 
study (i~8| . we include the in-medium ui (magenta-short- 
dashed line) and (light-brown-dash-dotted line) contri- 
butions to thermal di-electrons. The thermal di-electrons 
are dominated by the in-medium p mesons, while the u> 
contribution is submerged under the broadened p spec- 
tra. The thermal spectra with the CE EOS (green-dash- 
dotted-dotted line) has also been shown. The production 
rate in the PCE scenario is larger in the invariant mass 
range below free p mass than in the CE one, though the 
temperatures with PCE are lower. This is because the 
chemical potentials in the PCE scenario lead to a larger 
broadening of the p spectral function and an enhance- 
ment factor e 2 ^/ 7 " compared with the CE scenario. The 
enhancement at low masses in the PCE scenario is more 
obvious at lower Tf. We will come back to this issue 
later. These two EOS have the same partonic contribu- 
tions because their differences only occur in the hadronic 
phase. 

The dilepton emission rate from the freezeout vector 
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Figure 3: (Color online) (a) The invariant mass spectra of 
thermal di-electrons in full phase space. In the partonic 
phase, the main source is qq — > 7* — > e + e~ (magenta-long- 
dashed line). In the hadronic phase, the total contribution is 
the blue-dashed line, where the contribution from the p me- 
son (red-dotted line) dominates, and those from the uj and <f> 
mesons are shown in magenta-short-dashed and brown-dash- 
dotted lines respectively, (b) The invariant mass spectra of 
di-electrons for vector mesons p, uj and (j> at the freezeout. 



meson is given by 

dN a = «(tXl 

d*p 3 \ 9 J rv d 4 P ' [ ' 

where Ty is the total decay width of the vector mesons. 
The vector meson momentum spectra at thermal freeze- 
out can be expressed by the extended Cooper-Frye for- 
mula Q 
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Since the lifetimes of uj and <f> are much longer than the 
time scale of the freezeout process, we treat these con- 
tributions as in vacuum and neglect the medium effect. 
The imaginary parts of the uj and <fi propagators can be 
given by the Breit-Wigner formula, 
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Figure 4: (Color online) The cocktail mixture of the partonic, 
in-medium and freezeout hadronic sources for di-electrons. 
The partonic, in-medium and freezeout hadronic contribu- 
tions are in magenta-long-dashed, red-dotted, and blue-dash- 
dotted lines, respectively. The total contribution is in the 
black-solid line. 



But most of the p mesons decay in the medium due to 
its short lifetime, so we include the medium effect in the 
p meson propagator. In Fig.([3][b)) are shown the invari- 
ant mass spectra of di-electrons from the freezeout vector 
mesons, where the sharp peaks of the uj and <f> mesons 
can be seen compared to a much broader bump of the p 
meson. 

In Fig. we sum over all sources we have considered. 
The full mass spectra have two sharp peaks of the uj and 
<j) mesons at the freezeout same as in vacuum due to 
their long lifetime. For comparison the in-medium spec- 
tra are also shown where only a much lower peak from 
the <f> meson is visible, indicating clear medium effects. 
Subtracting these sharp peaks of uj and 0, the broad- 
ened spectrum of the in-medium p meson can be seen. 
The partonic contribution dominates over the hadronic 
one when M > 1.1 GeV/c 2 . These continuum-like IMR 
dileptons may provide a direct probe to the deconfine- 
ment phase transition in high energy heavy ion collisions. 
The different EOS give the similar structure but slightly 
different magnitude. It seems that the low mass enhance- 
ment favors the PCE scenario, we will come back to this 
point with details later in the following section. 



IV. COMPARISON WITH DATA 

Different from the SPS energy, the charm quarks have 
a considerable production rate at the RHIC energy. So 
there is a large background from semi-leptonic decays of 
the charm hadrons. In this section we will estimate this 
background and compare our dilepton results with the 
data. 

We use the event generator PYTHIA [34| (version 
6.416 with CTEQ5L PDF) to simulate the background 
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Figure 5: (Color online) The semi-leptonic decays of charm 
hadrons. (a) The re-scaled di-electron cross section from 
charm hadrons of semi-leptonic decays in p+p collisions by 
PYTHIA. The data are taken from the PHENIX collabora- 
tion (b) The nuclear modification factor for nonphotonic 
electrons in central Au+Au collisions from the PHENIX col- 
laboration [36j]. The fitting function is given in Eq.([8]). 



from semi-leptonic decays of the charm hadrons (D°,D , 
D s and A c ). The PHENIX collaboration also tuned the 
parameters of PYTHIA [35[ to fit the charm hadron data 
at SPS and FNAL and single electron data at ISR. The 
parameter dependences such as intrinsic kx and the par- 
ton distribution functions are also addressed in Ref. [14| . 
In our paper, we do not consider the fluctuations from 
these parameter. 

In p+p collisions, the dilepton yield in the mass range 
[1.1,2.5] GeV/c 2 is dominated by semi-leptonic decays of 
charm hadrons. In the PHENIX acceptance the inte- 
grated yield of di-electrons per event from heavy-flavor 
decays in that range is (4.21 ± 0.28 ± 1.02) x 10 -8 [11. 
With the branch ratio for charm quarks to electrons |37T | 
and the correction for the geometrical acceptance, the 
rapidity density of cc pairs can be estimated [l4j. We 
use the PYTHIA event generator with the PHENIX ac- 
ceptance to reproduce the spectra from charm hadron 
contribution, see Fig. ©a)). It can be seen that our 
results from PYTHIA (black-dashed line) are consistent 



with those given by PHENIX (red-dotted line). We ob- 
tain the cross section of cc pairs a cS = 0.5 mb. 

For Au+Au collisions, we use the renormalized cross 
section in pp collisions and scale it by the mean number of 
binary collisions. We choose N co u = 950 for most central 
collisions. The charm quarks are mostly generated in 
the pre-equilibrium stage. In medium the pt spectra of 
the charm quarks as well as the angular correlation of 
the cc pairs could be modified due to its interaction with 
the thermalized partrons. The medium modifications of 
heavy flavors have been widely studied in, e.g., Ref. [HI, 
[39( 1 . To include the medium modifications in a simple 
way, we parametrize the nuclear modification factor of 
the single electron in the form, 

R s aa(pt) = mux [1.0, exp (a/p T + b)} , (8) 

where a = 1.23 and b = —1.51, see Fig. JU^b)) for the 
fitting function and experimental data. Note that we 
have neglected the T contribution here. To get a realis- 
tic pt distribution for electrons, we use the original px 
spectra obtained by PYTHIA and multiply them with 
R aa (pt)- Using the Monte Carlo method, we sample 
the momentum spectra in accordance with the result- 
ing pt spectra for electrons and positrons respectively. 
In each event we randomly choose from the sample the 
momenta of one electron and one positron and combine 
them to a di-electron pair, from which the invariant mass 
and total transverse momentum of the pair can be deter- 
mined. The modified invariant mass spectra by R aa (pt) 
are found to be narrower than without such a modifica- 
tion as shown in Fig. (^a)). 

With the STAR acceptance (transverse momentum 
Pt > 0.2 GeV/c and pseudorapidity |?y e | < 1 for an indi- 
vidual electron, rapidity |y ee | < 1 for a pair of electrons) 
and pt resolution, we compute the di-electron spectra in 
most central collisions and compare them with the STAR 
preliminary data of 0-10% centrality, see Fig. ©a)). We 
also included the Dalitz decay channels for rj [40| and to: 
?y — > e + e~7 and u> — > e + e~ir° . The 77 contribution can be 
easily deducted as a background in the experiment due 
to its very long lifetime (about 1.5 x 10 5 fm/c), which 
leads to its decay outside the freezeout scope. 

The cooktail sum including the in-medium p mesons 
can roughly reproduce the di-electron spectra in the 
LMR, see Fig. @a)). In Fig. ©b)), we show the total 
di-electron spectra (thick lines), the contributions from 
the open charm (green-dashed line) and the in-medium p 
(thin lines) in the range M G [0, 1.2] GeV/c 2 . We found 
that the p meson contributions (thin black-solid line) are 
submerged under the open charm one. This indicates 
that the charm backgrounds play an important role in the 
dilepton spectra at the RHIC energy. As we discussed in 
Sect. |H| we can tune Tf to a lower value (e.g. 90 MeV 
in red-dashed line) to increase the contribution from the 
in-medium p mesons. The cocktail sum (thick red-dashed 
line) with the in-medium p contribution (thin red-dashed 
line) for Tf = 90 MeV is also shown. This seems to give 
a better fit to the data. Though the lower Tf gives larger 
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Figure 6: (Color online) The invariant mas s sp ectra and the 
comparison with STAR preliminary data [171 ] in most cen- 
tral (0-10%) Au+Au collisions with the STAR acceptance. 
The cocktail sums are in thick lines. See the text for de- 
tailed illustrations, (a) The results in M = [0.2, 3] GeV. (b) 
The results in M = [0.2, 1.2] GeV. The thin-balck/thin-red- 
long-dashed line denotes the contribution from the in-medium 
p decays for T f =106/90 MeV. The thin-red-dashed-dotted- 
dotted line denotes the contribution from the in-medium p 
decays without the pr cutoff for Tf =90 MeV. (c) The results 
in M — [1.1,3] GeV. The thin-balck/thin-red- long-dashed line 
denotes the QGP contribution for To =0.4/0.1 fm. 



broadenings of the p spectra and low mass enhancements, 
the p meson contribution is still smaller than the open 
charm one. This is because: (1) The nuclear modification 
factor enhances the charm contribution in the LMR; (2) 
Most of the low mass di-electrons from the in-medium 
p mesons have low px, which are beyond the capability 
of the detectors and can not be measured. To support 
the point (2), we calculate the in- medium p meson con- 
tribution incorporated by the STAR acceptance except 
the pt cutoff for electrons and positrons. The result is 
shown in the thin red-dashed-dotted-dotted line in Fig. 
@b)). We can see a strong enhancement below the free 
p mass. 

With the nuclear modification factor for charm 
hadrons, we can roughly reproduce the di-electron spec- 
tra in the IMR, see Fig. ((He)). One can see that the 
thermal contributions from the QGP phase (thin black- 
solid and red-dashed lines) are much smaller than the 
correlated charm decays (blue-dashed and green-dashed- 
dotted-dotted lines). Now we try to look at if it is pos- 
sible to increase the QGP thermal contributions in the 
IMR by the tuning parameters. We know that the ther- 
mal rate from the QGP is p roportional to T a At where 
A is the transverse area [T^. To this end, in our model, 
we can tune the equilibration time and entropy density 
(initial energy density) with the constraint sqTq = con- 
stant to keep the multiplicity rapidity density unchanged. 
The dilepton emission rates do not change much for dif- 
ferent ro. But for an earlier equilibration time, e.g. 
To = 0.1 fm (thin red-dashed line), the partonic con- 
tribution is enhanced in the IMR, since the early equili- 
bration time gives larger space-time volume of high tem- 
peratures, whose di-electron emissions mostly contribute 
to the IMR. But there is still a large gap between the 
contributions from charm hadrons and from the QGP. 
In addition, to lower the transition temperature will in- 
crease the space-time volume of the QGP phase and then 
dilepton rates from thermal partons. But this enhance- 
ment is almost in the LMR and will not significantly in- 
fluence the IMR. So it seems that it is very difficult to 
extract the thermal sources from the backgrounds from 
charm hadron decays in the invariant mass spectra alone. 
Additional observables such as pr spectra and collective 
flows fl8l . l4lT | are also needed. 

We show in Fig. © the results with the PHENIX ac- 
ceptance [ll| . We see that the charm backgrounds still 
out-perform the in-medium p. The acceptance geome- 
try pushes the charm hadron contributions toward the 
LMR. Using our cocktail sources, there is still a large un- 
expected excess of di-electrons in the LMR as reported 
by the PHENIX collaboration. 



V. SUMMARY AND CONCLUSION 

We investigate the di-electron low and intermediate 
mass spectra from the vector and charm hadrons in most 
central heavy ion collisions at ultra-relativistic energies. 
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Figure 7: (Color online) Comparison with PHENIX data Q]| 
in most central (0-10%) Au+Au collisions with the PHENIX 
acceptance. 

The space-time history of the fireball is provided by a 2+1 
dimension ideal hydrodynamic model, whose parameters 
are fixed by fitting the transverse momentum spectra of 
long-life hadrons, i.e., pions, kaons and protons. Two 
types of equations of state are used. The medium ef- 
fects of vector mesons from scatterings of vector mesons 
by mesons and baryons in the medium are considered. 
The di-electron emissions from in-medium vector meson 
decays can be evaluated via the imaginary parts of the 
vector meson propagators which are functions of space- 
time through the temperature. Due to their longer lives 
than the time scale of the freezeout process, most of the 
uj and mesons may decay at the thermal freezeout, giv- 
ing two sharp peaks in di-electron mass spectra. The 
contribution from the charm hadrons is modeled by the 
PYTHIA simulation of the proton-proton collisions and 



modified by the binary collision number and the nuclear 
modification factor for electrons. 

The cocktail sum over all above sources and the par- 
tonic phase incorporated with the acceptances of the 
STAR detector is compared to the STAR preliminary 
data. The hadronic many body effective theory with a 
broadening rho meson spectral function can describe the 
STAR di-electron data in the LMR. With a parametrized 
nuclear modification factor for electrons from charm 
hadron decays, we can roughly reproduce the di-electron 
spectra in the IMR, though we still lack enough knowl- 
edge about open charm decays in medium, such as the 
modification from the dynamical correlation of cc pairs. 

In conclusion, we find: (1) The detector acceptance 
especially the transverse momentum cutoff significantly 
suppresses the contribution from the in-medium p me- 
son in the mass region below the p mass; (2) With the 
current set of parameters and detector acceptances, the 
backgrounds from charm hadrons dominate in the low 
and intermediate mass regions. Therefore it is impos- 
sible to extract the thermal sources of dileptons with 
the invariant mass spectra alone if the backgrounds from 
charm hadrons are not removed. Other observables such 
as transverse momenta and collective flows may provide 
additional tools to tag these sources. Future STAR pro- 
grams such as the Heavy Flavor Tracker [12] and the 
Muon Telescope Detecctor are expected to improve the 
capability of identifying the backgrounds from charm de- 
cays and extracting the thermal souces. 
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